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1. Introduction 



The discussion of electronic correlation effects in 3rf-transition-metal compounds appears 
often dominantly reserved for the physics of oxides, e.g. for cuprates, manganites, 
5^ \ ruthenates, or more recently cobaltates. However there is increasing awareness 

that strong correlation effects are also relevant for many intriguing electronic phase 
transitions of various chalcogenides such as the metal-insulator transition (MIT) 
in NiS2-xSea; [Ij, the charge- density wave (CDW) transition in NbSe3 [2j or the 
recent discovery of higher-temperature superconductivity in Fe- chalcogenides. In this 
respect, the vanadium chalcogenides with the quasi- two-dimensional X-VS2 misfit-layer 
structures [3j and the quasi-one-dimensional (ID) BaV(S,Se)3 systems offer a unique 
opportunity to study the interplay of correlation effects and effective low- dimensionality 
in these transition-metal compound systems. 

Here we want to concentrate on the subtle comparison of the quasi- ID barium 
vanadium chalcogenide systems. It will become evident that the vanadium selenide 
BaVSea is in effect rather similar to the corresponding vanadium sulfide BaVSa iii many 
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respects. The latter compound has been the subject of various studies concerning its 
complex electronic structure throughout the last decades [4l-[T2]. This sulfide exhibits 
three continuous phase transitions with decreasing temperature, one of which is a 
highlighting metal-to-paramagnetic-insulator transition at Tmit ^ 70K. Though the 
exact details of the scenario are still to be revealed, strong evidence exists [IIl[l3] that 
this transition is associated with a CDW instability. Interestingly, such a transition 
cannot be found in the similar selenide compound, which remains metallic down to 
very low temperatures [14j. This raises questions about the delicate electronic and 
structural differences that anticipate a possible CDW in these chalcogenides. Hence 
the present investigation is devoted to a detailed comparison between the correlated 
electronic structure of BaVSa and BaVSe3. 

Concerning the structural properties, many features of BaVSea are akin to the 
ones of BaVSa. At room temperature, both materials exhibit a hexagonal structure 
in which the vanadium atoms form straight chains in the crystallographic c direction. 
The intrachain V-V distance is less than half of the interchain distance, so that one is 
tempted to identify a quasi- ID substructure in these compounds. The chains become 
slightly distorted in the ab plane with a zigzag distortion in the be plane at the 
continuous transition temperature Ts ^ (240K, 290K) (sulfide, selenide) p!5lfT6] . 
Thus, below this temperature regime, an orthorhombic crystal structure is stabilized for 
both materials. Furthermore, the additional metal-to-insulator transition of the sulfide 
is accompanied by a doubling of the unit cell and a transition to the monoclinic crystal 
system [13j, whereas the selenide remains orthorhombic down to very low temperature. 
Both chalcogenides display local-moment behavior at ambient T and a final continuous 
magnetic ordering transition. While in BaVSa ctn incommensurate antiferromagnetic 
(Neel) order develops below Tn ^ 30K [TTj, BaVScs undergoes a ferromagnetic 
transition at Tc 43K fTGlfTS]. 

The electronic structure of both materials in the low-energy regime is dominated by 
the nominal valence, so that they can be described as 3d^ compounds with a t2g-like 
manifold at the Fermi level |15[[l0l[12j . The latter consists of an yli^-like orbital pointing 
along the vanadium chain direction as well as two £'^-like orbitals pointing inbetween 
the neighboring sulfur /selenium ions, respectively. Due to strong hybridization with 
S(3]9)/Se(4]9), the V(e^) states have mainly high-energy weight and may be integrated 
out in a first low-energy approximation. Hence the essential physics of these systems 
can be described within a three-band model approach [12j. The two Eg states are 
degenerate in the hexagonal structure, but split into an Egi and an Eg2 state in the 
orthorhombic regime. Due to their somewhat isolated orbital space, the hybridization 
with the environment is comparably small. For this reason, they form very narrow bands 
that show almost no dispersion. In contrast, the Aig orbitals display a larger dispersion 
along the V chain direction. Hence, in a minimal model, the physics of these materials 
essentially boils down to the interplay of an itinerant and two localized states sharing a 
single electron. 
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2. Theoretical Framework 

Due to the obviously delicate electronic systems in BaVSa and BaVSea, with narrow 
bands at the Fermi level and local-moment physics, methods are needed that allow 
us to take into account many-body effects explicitly, i.e. beyond the standard 
density functional theory (DFT) in its local density approximation (LDA) to the 
exchange-correlation energy. In the present work, the opportunity is taken to compare 
two different approaches for this task, namely the LDA+DMFT p^[20j framework 
in conjunction with a highly- evolved quantum-Monte-Carlo impurity solver as well 
as the combination of LDA with the rotationally-invariant slave-boson formalism 



The DFT parts of the following calculations have been performed using 
a mixed-basis pseudopotential (MBPP) code [23], which utilizes norm-conserving 
pseudopotentials and a combined basis set of plane waves as well as non-overlapping 
localized functions. From the LDA band structure, by means of the maximally-localized 
Wannier function construction j24l - [26] . a low-energy Kohn-Sham (KS) Hamiltonian 
i7^g(k) within a correlated subspace C is extracted [27j. In this work we restrict all 
self-energy effects to this subspace and hence may write the local interacting Green's 
function as 



Here ujn are Matsubara frequencies and S^^^ is the self-energy matrix for the correlated 
orbitals denoted by m^m' . Note that ([T]) is written in the DMFT approximation, i.e., 
with a purely local self-energy (for a review see e.g. [28j). The mean- field version of 
RISB which is put into practice in the present work also has a local self-energy by 
construction. Moreover, the RISB self-energy only captures the linear term of the 
frequency dependence, expressed via the quasiparticle weight Z, as well as a local 
static part [12j. Thus without going into details, the main differences between DMFT 
and RISB (in the mean-field version) is the neglection of the frequency dependence 
of the higher-energy excitations in the latter framework. However, as a saddle-point 
approximation the performance is fairly efficient, with reliable qualitative (and often 
even good quantitative) results in most cases. It is also important to realize that 
although we use here a lattice-implementation of the RISB method, this technique may 
equally well be utilized as an impurity solver for standard DMFT. Due to the inherent 
local nature at the saddle-point, this specific DMFT impurity solution is then identical 
to the direct lattice-calculation result. 

In the present work, the explicit quantum impurity problem within DMFT is 
solved using the continuous-time quantum-Monte-Carlo (CTQMC) approach employing 
the hybridization-expansion method [29]. For the actual computations the recent 
implementation by Ferrero and ParcoUet (see [SQlEI]) is employed. Both frameworks, 
i.e., DMFT (CTQMC) and RISB, allow for non-density-density terms in the many- 
body Hamiltonian, which we restrict in the present case to spin-flip and pair-hopping 
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processes. Thus the following interacting Hamiltonian for the minimal modeling may 
be used: 

Hint = UY^ flm^flmi + ^ flmcjflm'a + {U' - J) flmcjnm' a] 

where n^cr = ^ma^ma^ with rf^)^ as the electron annihilation (creation) operator for 
orbital m and spin g. Throughout the calculations, the choice U' = U-2J and Jq = J-> 
appropriate for t2g systems, is made. 



3. Results 

3A. LDA studies 

The following LDA comparison between BaVSa and BaVSea builds up on the 
orthorhombic crystal structure (space group (7mc2i), since we are mainly interested 
in elucidating the possibility of the CDW state in both systems. Note that the CDW 
instability is not only absent in the selenide compound, but also vanishes in the sulfide 
at high pressure [9j. This raised the suscpicion that the selenide compound may be 
interpreted as the high-pressure analog of the sulfide compound [14j. In any case, 
theoretically investigating BaVSea should surely also provide more information about 
the origin of the apparent CDW instability in BaVSa. 

The existence of the orthorhombic phase of the selenide has been reported already 
on the occasion of its first synthesis by Kelber et al. [16j, but to the present authors' 
knowledge no detailed experimental structural data is available in the literature. For 
this reason, the atomic positions used for (7mc2i-BaVSe3 are revealed by structural 
relaxation within the LDA using the MBPP code, while the crystal parameters are based 
on similar recent calculations performed by Akrap et al. [14j. In contrast, structural 
data obtained by neutron diflFraction by Ghedira et al [l32j are used for BaVSa. Note 
that an LDA structural relaxation of this experimental data for (7mc2i-BaVS3 does 
result in only marginal changes of the atomic positions, therefore an influence thereof 
on the following studies cannot be found. The crystal parameters (a, 6, c) for one unit 
cell consisting of two formula units of BaVSa or BaVSe3, respectively, in au read (12.77, 
2L71, 10.58) for BaVSa and (13.34, 22.68, 11.06) for BaVSea. For both materials 
the ideal quasi- ID chains of vanadium ions are slightly distorted in the orthorhombic 
phase, leading to zig-zag chains. This effect turns out to be slightly larger in the 
selenide (relative displacement: 0.025) than in the sulfide (0.021). Concerning the 
electronic structure, figure [H shows the LDA density of states (DOS) of BaVSes in a 
comparatively large energy window around the Fermi level. Analogous to BaVSa [E], 
it can be seen that the low-energy physics of the material may be dominantly described 
by the abovementioned t2g manifold. Thus a downfolding to an effective three-band 
model consisting of a broader 74i^-like level and narrower E'^i-like and £'^2-like levels is 
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Figure 1. LDA density of states of BaVSes. 
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Figure 2. Left: LDA band structure of BaVSs (top) and BaVSes (bottom). The 
fatband broadening depicts the weight of t2g orbitals on the bands. Color coding: 
Aig^ blue/dark, Egi^ red/medium and Eg2^ green/hght. Right: Brihouin zone of the 
orthorhombic structure with relevant points/directions. 

a suitable approximation for the essential physics. The accordance with the LDA DOS 
of the sulfide (shown in [12]) on this level of comparison is rather striking and only minor 
differences show up. Figure [2] compares the band structure of BaVSca and of BaVSa in 
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a smaller energy window, with 'fatbands' showing the weight of the symmetry-adapted 
t2g states on the specific bands. Again both band structures are similar, yet subtle 
diflFerence may be observed. First the Se(4]9) band at the F point is much closer to the 
Fermi level than the corresponding S(3p) band in BaVSa, which may lead to a smaller 
Hubbard U for the minimal 3-band model due to increased screening. Furthermore 
the bandwidth difference between the three t2g bands is reduced, i.e., the yli^-like bands 
become narrower and the £'^-like ones broader. Hence the total t2g bandwidth in BaVSe3 
is slightly smaller (^2.4 eV) than in BaVSa (^2.7 eV). This may be explained by the 
fact that since the larger Se ions give rise to increased interatomic distances, a smaller 
74i^-like hopping amplitude along the chains results. As a general feature, the overall Aig 
contribution to the bands seems somewhat more distributed in BaVSca. This enhanced 
spread in energy is again providing hints towards a larger degree of, eventually more 
isotropic, localization in the selenide. 




Figure 3. LDA Fermi surfaces of BaVSs (top) and BaVSes (bottom) from different 
perspectives. 

Experimentally the nesting condition for the formation of the CDW was detected 
as 2kF = 0.5c* [H], where c* is the reciprocal lattice vector along the c axis. In order 
to theoretically examine this condition, it is necessary to study the topology of the 
Fermi surface (FS) (see figure [3j). Therefrom the nesting condition is not fulfilled in 
the pure LDA result in both materials. However somewhat surprisingly, the quasi- 
ID sheet of the selenide is closer to fulfilling the experimental nesting condition than 
the corresponding sulfide sheet in LDA. In contrast, the pillar-like structures are much 
more enhanced in the selenide. Hence differences in the t2g occupations between the two 
compounds may be expected. Figure [4] shows the respective three-band dispersion of 
the correlated subspace obtained from a maximally-localized Wannier construction [24j 
for strongly hybridized bands [25] . It should be mentioned that this construction yields 
for the selenide a larger spread of the Wannier functions (i.e. in the range of 23 au^ 
for the £'^-like ones in BaVSe3 compared to about 18 au^ for the same ones in BaVSs). 
This can be attributed to a stronger hybridization of the t2g system with the Se(4]9) 
electrons. Note that for the further use of this Wannier representation we rotated the 
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Figure 4. Dispersions of the effective three-band model (red/gray), obtained by 
means of a maximahy-locahzed Wannier-function construction, for BaVSs (top) and 
BaVSes (bottom). Black: original LDA band structure. 

original maximally-localized Wannier Hamiltonian into the crystal- field basis, i.e., the 
basis where the on-site part thereof takes on a diagonal form. 

Concerning the susceptibility towards the CDW state it is important and evident 
to note that the nesting vector mentioned above is not necessarily located along the T-Z 
high-symmetry line jl2il27j. In order to take this fact into account, the low-energy band 
structure is computed along a Brillouin zone (BZ) path containing a line parallel to the 
T-Z line, shown in figure [5l The points 'M/2' and 'A/2' are located halfway between 
the F- or Z-point, respectively, and the edge of the first BZ (see figure [2j). From the 
displayed fatbands it can be seen that the hybridization between Aig and Egi on the 
resulting bands along M/2- A/2 as well as along V-Z is significantly larger in the selenide 
than in the sulfide. Thus a clear notion of exclusive ''Aig bands' or ''Egi bands' becomes 
even more difficult for BaVSca. This fact also impedes the assignment of the distinct 
FS sheets to a single definite orbital character. 

3.2. Inclusion of correlation effects 

Besides the apparent metal-to-insulator transition in BaVSs, already the complex 
magnetic behavior in both systems (with different magnetic order at low temperature) 
shows clearly the importance of electronic correlations in these vanadium chalcogenides. 
Hence a deeper comparison of BaVSs and BaVScs has to include explicit many-body 
physics, which we approached with the LDA+DMFT (CTQMC) method as well as 
the LDA+RISB technique. This twofold investigation was elaborated in order to 
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Figure 5. LDA t2g Wannier fatband plots for BaVSs (top) and BaVSes (bottom) 
along a closed path in the first Brillouin zone, connecting F, 'M/2', '^4/2' and Z (see 
figure [2]). Color coding: Aig (blue), Egi (red) and Eg2 (green). 



evaluate the reliability of the latter method compared to the more complete CTQMC 




Figure 6. (a), (b): Effective t2g fillings for BaVSs (dashed/squares) and BaVSes 
(solid/circles) using the KS- Wannier Hamiltonians with increasing U. (c), (d): Fillings 
for BaVSes using only density-density terms (dash-dotted/triangles) and including 
spin- flip as well as pair- hopping terms (solid/circles), (a), (c): CTQMC solution of 
DMFT and (b), (d): MSB solution. 
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impurity solution within DMFT. In the foUowing, aU the presented CTQMC impurity 
computations were performed at T = 116K (/3 = 100 eV~^) and including spin-flip and 
pair-hopping terms in the interacting hamiltionian (if not stated otherwise). 

In [33j, a chosen value of [/ = 3.5 eV led to a proper theoretical description of 
BaVSa. Iii this work we also do not compute the interaction parameters from flrst 
principles, but again take a practical point of view in order to account for the key 
correlation eflFects. As noted in the previous section, the corresponding Hubbard U 
for BaVSe3 is expected to be smaller because of the increased screening capabilities 
of the Se(4]9) electrons compared to the S{3p) ones. In order to compare BaVSa and 
BaVSe3 in this respect, the calculated orbital-resolved flllings are shown in flgure[6]for 
varying values of U and flxed Hund's exchange J = 0.7 eV. First note the very good 
agreement between the DMFT (CTQMC) and the RISB treatment of the correlation 
eflFects. Both methods clearly show the correlation-induced charge transfer mainly 
between the Aig and Egi orbitals [I2[[33]. The actual numbers surely differ somewhat 
between both approaches due to the neglect of the explicit quantum fluctuations within 
RISB, especially concerning the frequency dependence of the high-energy excitations. 
But also the trend in the (minor) diflFerences between density-density only terms and the 
inclusion of more general terms in the interacting Hamiltonian is still well reproduced by 
the simplifled RISB technique. Concerning the physics, it is visible that the Aig fllling 
is generally reduced in the selenide compared to the sulflde, which was already expected 
from the Fermi-surface discussion. This result is also apparent from table [U where the 
occupation numbers are displayed for selected interaction strengths. In this context, 
a smaller Hubbard interaction U = 2.5 eV for BaVSea seems reasonable in the end. 
However, for direct comparisons we also include the case U = 3.5 eV in the computations. 
Note that J is usually not strongly affected by the crystal environment and to a good 
approximation can remain constant. Hence, if not stated differently, we generally flxed 
the Hund's coupling to the value J = 0.7 eV used in previous studies [12j and which 
is in line with what was utilized in similar investigations for vanadium oxides [341436] . 
Note that a suitable energetical range for J proved to be sufficient to describe the key 
low-energy physics of BaVSa [33j . 

In summary, the system could be understood as an effective two-band system, since 
the effective Eg2 orbital is almost empty, with an even larger correlation-induced charge 
transfer from Aig to Egi in the selenide than in the sulfide. It can further be seen that the 
spin-fiip and pair-hopping terms have only a small infiuence on the occupation numbers, 
resulting qualitatively in a slightly smaller orbital polarization. The latter is due to the 
fact that those terms lead to a further reduction of the integrated effective U compared 
to the density-density only case. Moreover since the present systems mainly reside in the 
local single-particle sector, the infiuence of spin-fiip and pair-hopping matrix elements 
on the charge fiuctuations remains minor. 

Figure [71 shows the density of states derived from the original Wannier functions 
compared to the local spectral functions from DMFT for both systems. On the LDA 
level, the plot shows the reduced bandwidth of the selenide compound, with the orbital 
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Table 1. Orbital-resolved fillings from LDA+DMFT (CTQMC). Vanishing Hubbard 
U marks the LDA result. 'DD': only density-density interactions, 'RF: additional 
inclusion of spin-flip and pair-hopping interactions. 





U (eV) 


^int 






Eg2 


BaVSs 


0.0 




0.58 


0.31 


0.11 




3.5 


DD 


0.38 


0.54 


0.08 




3.5 


RI 


0.37 


0.53 


0.10 


BaVSea 


0.0 




0.49 


0.40 


0.11 




2.5 


DD 


0.40 


0.46 


0.14 




2.5 


RI 


0.41 


0.45 


0.14 




3.5 


DD 


0.31 


0.63 


0.06 




3.5 


RI 


0.33 


0.61 


0.06 



.-2-10 12 3 



1 1 2 3 4^ 




Figure 7. Comparision of the LDA Wannier DOS (left) and local spectral functions 
from LDA+DMFT (right) for BaVSs with U = 3.5 eV (top) and BaVSes with U = 2.5 
eV (bottom). 



resolved effective widths within the t2g manifold becoming more similar. Furthermore, 
the overall stronger hybridization (note the use of the crystal-field basis) between 
Aig and Egi in BaVSes is clearly visible. This also leads to a more pronounced Aig 
quasiparticle peak right at the Fermi level in the DMFT local spectral function. 

As explained in [12l|33] , the importance of explicit many-body eflFects in deforming 
the LDA Fermi surface appears as a crucial ingredient for the onset of the CDW in 
BaVSs- Therefore, figure [8] shows the DMFT (CTQMC) quasiparticle Fermi surfaces 
for both chalcogenides. Correspondingly, figure [9] displays the quasiparticle dispersion 
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Figure 8. Quasiparticle Fermi surfaces of BaVSa with [7=3.5 eV (top) and BaVSes 
with [7=3.5 eV (middle) and [7=2.5 eV (bottom). 

for the closed path in the first BZ connecting the points 'M/2' and '^4/2' (see figure [2j), 
which mark the relevant distance for the possible CDW nesting. The plots render it 
obvious that no definite information about the absence of the CDW in BaVSea can be 
retrieved from this viewpoint. Since qualitatively the same charge-transfer mechanism 
applies to the selenide compound, it may be observed that the FS deformations also 
go in the same direction. This means that within the LDA+DMFT accuracy also for 
smaller Hubbard U the experimental nesting condition could in principle be realized in 
BaVSe3. Note however that in the selenide for [7=2.5 eV the respective quasi-lD sheets 
are still not strongly flatted. 

In order to get a quantitative comparative understanding of the degree of quasi-one- 
dimensionality of the two materials, we investigated the theoretical order of magnitude of 
the anisotropy in the DC conductivities in LDA as well as in the strongly correlated case. 
To this end, table [2] summarizes the relation of the spatial components of the averaged 
Fermi velocities squared ((v|.) = {vD + {Vy) + ('^^)), whereas y denotes the spatial 
direction in which the zig-zag distortion of the vanadium chains can be observed. Note 
that the LDA values obtained for the sulfide are in perfect agreement with the earlier 
result from Mattheiss [37] . Surprisingly, the selenide shows an even larger formal DC 
anisotropy in LDA compared to the sulfide compound. In contrast, it is also obvious that 
within the many-body approach of LDA+DMFT the renormalization of the anisotropy, 
especially in the dominant {vD / {v'^) channel leads to a balanced value for BaVSe3 with 
reasonable U=2.5 eV. In contrast, the quasiparticles become indeed significantly heavier 
in the c direction for BaVSa with U = 3.5 eV. It is therefore worthwhile and interesting 
to note that the sulfide compound indeed displays a clear signature of DC anisotropy in 
the correlated case, whereas such a clear character is appearing in the according BaVSe3 
data only at improper large U. It shall be remarked that the intraorbital quasiparticle 
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Figure 9. Effective t^g Wannier band structure with (red/gray) and without 
(black/dark) renormahzation from LDA+DMFT along the closed path including 'M/2' 
and 'A/2' (see figure [2j) for BaVSa (top) and BaVSes with /7=3.5 eV (middle) and 
t/=2.5 eV (bottom). 



Table 2. Relation of the spatial components of squared Fermi velocities in BaVSs and 
BaVSea on the LDA level and within LDA+DMFT (CTQMC) for a given V (in eV). 



Ratio 



BaVSs 
LDA [7=3.5 



3.7 
8.3 



2.7 
9.9 



LDA 

4.8 
10.7 



BaVSea 
[7=2.5 U-- 

5.2 
9.9 



:3.5 

3.0 
10.2 



weight for BaVSs (BaVSea) amounts to Z ^ 0.5 (0.6) with minor differences within 
the t2g manifold. Note that, although we work in the crystal-field basis with a diagonal 
onsite KS Hamiltonian for each of the two V atoms in the unit cell, the interacting 
terms together with the hybridization between these atoms lead to additional offdiagonal 
(interorbital) self-energy terms between Aig and Egi for the present crystal symmetry. 
Both many-body approaches, i.e. DMFT (CTQMC) and RISB, are capable of revealing 
and handling this eflFect. Thus a diagonal Green's function approach is invalid in the 
present case. 

Finally, since the spin degree of freedom plays an additional vital role in these 
chalcogenide compounds, figure [lO exhibits the respective on-site spin-correlation 
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Figure 10. On-site spin correlations for BaVSs (left) and BaVSes (right). Dashed 
lines: density-density only interactions, solid lines: including also spin- flip and pair- 
hopping terms. In each U scan a ratio U/J=b was chosen, (a-b): Total and orbital- 
resolved diagonal spin-spin expectation value, (c-d): Off diagonal spin-spin correlation 
between the effective t2g orbit als. 

functions from the LDA+RISB computations for constant ratio U/J=5. Concerning 
the diagonal spin-spin expectation value, as expected, the magnitudes evolve according 
to the nominal orbital occupations. Hence {S'^) for the effective Egi orbital dominates 
over the corresponding value for the effective Aig orbital in the strong-interaction 
limit. In this respect, again quantitatively BaVSes exceeds BaVSs due to the increased 
Egi occupation. The total {S'^) increases with larger U due to the growing electron 
localization. Note that in the present problem the magnitude of this diagonal quantity 
can exceed the atomic-limit value of a single electron, i.e., (S'^)at=|(| + = albeit 
we have nominally one electron in the 3d shell of the chalcogenides. This is due to 
the fact that at some point a strong Hund's coupling may tolerate an occupation of 
two electrons in the local three-orbital manifold of a metal. Since two electrons with 
parallel spin alignment form an 5^=1 system, the resulting {S'^) is bigger than the sum 
of two individual single-electron occupied sites. The interorbital spin-spin correlations 
again prove the dominant Aig-Egi hybridization, which also leads to a significant spin 
coupling. However, it may be observed that in the minimal modeling for both systems 
the first Hund's rule is violated close to the [/=0, i.e., the offdiagonal {SAigSsgi) becomes 
smaller than zero in this regime. 

It also has to be noted that these interorbital spin-spin correlations are rather 
sensitive to the degree of approximation of the rotational invariance in the interaction. 
While the diagonal spin-spin expectation value is not strongly affected by the inclusion 
of spin- flip and pair-hopping terms in the local Hamiltonian, the offdiagonal {SmSm') is 
strongly enhanced in this case by a factor of two to three. Thus, whereas integrated one- 
particle quantities like occupation numbers and local spins do not strongly suffer from 
density-density only descriptions, the two-particle functions, i.e. susceptibilities, may be 
substantially different. Figure [TT] displays the same spin-spin functions, but for the (more 
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Figure 11. Same as figure [TOl but with fixed Hund's coupling J =0.7 eV. The vertical 
dotted lines indicate the expected suitable choice for the Hubbard U of the respective 
system. 



realistic) case of a constant J=0.7 eV. Of course, the strongest response may be observed 
for smaller U since then the Hund's coupling is most effective. Note again the differences 
for the density- density only interactions, leading to an even qualitatively different 
behavior for small [/, where {SAigSEgi) is eventually dominated by the {SEg2SEgi)- In 
the more sound close to rotationally invariant description however, {SAigSsgi) is always 
clearly strongest. Because of the supposedly different U/ J ratio for BaVSa and BaVSca 
we hence expect for the selenide compound a more distinct magnetic behavior. Note 
that the interorbital spin-spin correlations decay in any case for large [/, contrary to 
what is known for half- filled systems (there {SmSm') actually increases with U). This is 
due to the fact that the ground state for a multi-orbital problem with a single-electron 
filling corresponds at large U to the localization of that electron within a single orbital. 

4. Discussion 

From the above calculations, it turns out that BaVSa and BaVSea are in many respects 
very similar materials. The differences between them are rather subtle, but still 
large enough to be resolved by the applied methods. The picture that arises thereof 
is not totally unambiguous, but several hints can be given why the formation of a 
charge-density wave is hindered in BaVSea. Already on an LDA level, it can be seen 
that hybridization effects among the bands of the t2g multiplet as well as with the 
corresponding S{3p) /Se{4:p) electrons, are significantly stronger in the selenide than in 
the sulfide. Hence a manifest discrimination between dominant Aig and dominant Egi 
character is far less obvious than in BaVSa. Furthermore, the strong role of the effective 
Aig orbital is weakened already on the LDA level and the final dominance of the effective 
Egi orbital in the correlated case is even enhanced. This is altogether underlined by the 
DC anisotropics for the selenide, that under the effect of the expected correlations do 
not show the highlighting renormalizations in the c direction. Thus although the orbital 
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differences are obvious, it seems as if the t2g manifold in BaVSea acts somewhat more 
cooperatively compared to the largely competing scenario in BaVSa. 

With a surely smaller Hubbard U the selenide compound should also be closer to 
the LDA limit [Hj, with a less pronounced nesting susceptibility towards the CDW 
state. However as described, from a pure Fermi-surface discussion no clear answer to 
this question may be conveyed. Nevertheless, electronic correlations are not irrelevant 
in BaVSe3 as they underline the dominance of the Egi level and the eventual importance 
of magnetic correlations. These latter are again stronger in the selenide compared to 
the sulfide, at least from a local viewpoint. Further calculations of intersite spin-spin 
correlations should shed more light onto this. 

In the end, the present study within standard LDA+DMFT is not in the definite 
position to clearly assign the diflFerent experimental facts of these suprisingly similar 
compounds from a theoretical perspective. However, although the situation is very 
subtle, we believe we have shown several indications that may motivate to a certain 
extent why the present compounds behave so diflFerently at lower temperatures. There 
seem to exist two main possible instabilities in these chalcogenides, namely the CDW 
one and the (ferro) magnetic ordering instability. As the driving forces towards these 
broken-symmetry states compete, a tailored renormalization-group study could reveal 
the finally dominating instability channel for each compound. Hence it is expected 
that the minor differences revealed in the present investigation show up more evidently 
within a renormalized scaling scenario. Still it has to be appreciated that, though not 
fully comprehensive, the LDA+DMFT method is nowadays on a level where it is feasible 
to reveal fine details between rather similar materials. 
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